Controlling Energy Gap of Bilayer Graphene by Strain 
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Using the first principles calculations, we show that mechanically tunable electronic energy gap 
is realizable in bilayer graphene if different homogeneous strains are applied to the two layers. It is 
shown that the size of energy gap can be simply controlled by adjusting the strength and direction of 
these strains. We also show that the effect originates from the occurrence of strain-induced pseudo- 
scalar potentials in graphene. When homogeneous strains with different strengths are applied to each 
layer of bilayer graphene, transverse electric fields across the two layers can be generated without 
any external electronic sources, thereby opening an energy gap. The results demonstrate a simple 
mechanical method of realizing pseudo-electromagnetism in graphene and suggest a maneuverable 
approach to fabrication of electromechanical devices based on bilayer graphene. 

PACS numbers: 73.22.Pr, 77.65.-j, 62.25.-g, 81.05.ue 



Materials that can change energy gaps under mechani- 
cal deformations are greatly desired to develop electrome- 
chanical nanodevices with controllable optical and elec- 
trical operational ranges [llM- Graphene, one atom thick 
semi- metallic membrane [3-5], may have this property 
because of its special electromechanical characteristics 
such as strain-induced phases in the chiral wavefunctions 
of the massless Dirac Hamiltonian (El-Hoj. Although ex- 
ternal electric fields have been used to tune the energy 
gaps in several forms of graphene ll-2o|, a simple and 
practical method of tuning the electronic energy gaps 
based on mechanical operations is lacking. 

The influence of smooth elastic deformations on the 
physics of a single layer of graphene can be described by 
introducing a suitable gauge-field vectorpotential to the 
free massless particle Dirac's equation [5j. The choice of 
vector potential depends on the direction and strength 
of applied strains. Specifically designed inhomogeneous 
strains can produce a net pseudo-magnetic flux that be 
expected to induce interesting electromagnetic effects [6|- 
llQI |. In contrast, under homogeneous uniaxial or biaxial 
strains, the effect of pseudo-magnetic fields is absent but 
a homogeneous pseudo-scalar potential can be generated 
owing to electron density variations or dilations EEL 
I23I ]. Because graphene is a one atom thick semimetal, the 
pseudo-scalar potential will alter its work function signif- 
icantly. When a uniaxial (or isotropic) strain is applied 
up to an experimentally realizable magnitude of ~10%, 
the work function of graphene is predicted to increase by 
0.27 (or 0.65) eV from the value of 4.49 eV of graphene 
in equilibrium [23|. Hence, the homogeneous strain can 
have important applications in various electronic devices 
for tuning the band lineup between metallic leads and 
graphene, or in changing the amount of charge trans- 
fer from foreign molecules. We also expect that such 
strain-induced pseudo-scalar potentials will modify the 
electronic structure of bilayer graphene significantly. 



In this paper, we predict the formation of mechani- 
cally tunable electronic energy gap in bilayer graphene. 
It is shown that the perpendicular electric fields across 
the two layers of bilayer graphene can be generated with- 
out any external gate potential if each layer is subjected 
to different strengths of homogeneous strains. This ef- 
fect originates from asymmetric generations of pseudo- 
scalar potential in each layer of strained bilayer graphene, 
which alter their work functions significantly [13]. We 
also show that the strain-induced electronic energy gaps 
are indirect in general and strongly depend on the size 
and direction of applied homogeneous strains. We dis- 
cuss origins of such peculiar variations of the bandgaps 
and present realistic experimental setups for observing 
the energy gaps. 

Our study of electronic structures of strained bilayer 
graphene is based on the pseudopotential density func- 
tional method [24J. The exchange-correlation interac- 
tions were treated within the local density approxima- 
tion [25]. The cutoff energy for expansion of wave- 
functions and potentials was 400 eV. Extensive momen- 
tum space sampling in the whole Brillouin zone (BZ) 
is performed to obtain accurate energy dispersions, re- 
liable energy gaps, and their dependences on strain. The 
Monkhorst-Pack /c-point grid of 2 x 12 x 1 is used for 
the atomic relaxation and of 5 x 30 x 1 for electronic 
structure calculation of 11.1% asymmetrically strain bi- 
layer graphene. Calculations for other strengths of strain 
also used an equivalent number of /c-points. The atomic 
relaxation was carried out until the Helmann-Feynman 
forces were less than 0.03 eV/A. The distance between 
adjacent bilayer graphene in the supercells is 26 A avoid- 
ing spurious interactions between them. We also have 
increased the distance up to 36A finding no changes in 
electronic structures. The weak dispersion force between 
the two layers of bilayer graphene is not considered here 
and the interlayer distance is fixed to 3. 35 A because its 
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FIG. 1: (Color online) (a) Bilayer and single layer graphene 
junction between the source and drain. Red arrows: ten- 
sile force direction, (b) Schematic atomic structure of bilayer 
graphene at equilibrium, showing carbon atoms in the top 
layer (clear circles) and the bottom layer (green circles), and 
the bilayer graphene unit cells with four carbon atoms (grey 
parallelogram) and eight carbon atoms (blue rectangle), (c) 
Atomic structure of asymmetrically strained bilayer along the 
zigzag chain direction (AZS) with 11.1% strain on the bottom 
layer, (d) Asymmetrically strained bilayer in armchair direc- 
tion (AAS) with 20% strain on the bottom layer, (e) Top 
panel: the first Brillouin zones (BZ) of bilayer graphene in 
equilibrium for the unit cell with four carbon atoms (black 
hexagon) and for the unit cell with eight carbon atoms (blue 
dotted rectangle). Bottom panel: first BZ for the AZS in (c) 
(red dotted rectangle); first BZ for AAS bilayer in (d) (green 
dotted rectangle). 



effect is negligible when obtaining the energy bands using 
the present calculation scheme [26| . We note that present 
methods correctly describe the trend of the energy gap 
in bilayer graphene under external electric fields with a 
typical underestimate of the gap size compared with a 
recent optical measurement [20]. 

Homogeneous strains can be applied to the two lay- 
ers of bilayer graphene asymmetrically in various device 
configurations. The prototypical system considered here 
is a junction between single-layer graphene and bilayer 
graphene [Fig.QJa)]; this junction can be synthesized us- 
ing the graphene transfer technique [I?}, and many me- 
chanically exfoliated samples contain such junction ge- 
ometries naturally {is). If a uniform tensile force is ap- 
plied to the junction [Fig. QJa)], each layer in bilayer 
graphene experiences a different magnitude of homoge- 
neous strain, because only one end of the top layer in 
bilayer graphene is pinched by the metallic lead whereas 
the other end is free-standing [Fig. Ufa)]. In this case, the 



top layer of bilayer graphene in the junction maintains its 
equilibrium atomic structure due to the negligible inter- 
layer interaction along the in-plane direction [29| , but the 
bottom layer is subjected to biaxial strain because both 
of its ends are clamped to the stretched substrates by the 
metallic leads [Fig. Ufa)]. Generally speaking, any device 
structure with asymmetric contacts to the top and bot- 
tom layers can produce asymmetric strains in each layer 
due to external tensile forces. Hereafter, we apply the pe- 
riodic boundary condition to the asymmetrically strained 
bilayer graphene only. To satisfy the commensurability 
requirement for between the top and bottom layers along 
the direction of applied tensile force, the unit cell of the 
system must be a rectangle with a very high aspect ratio 
[Fig. [He) and Hfd)]. Accordingly, the first BZ of the 
strained graphene is folded into a small rectangle within 
the first BZ of the original bilayer graphene in equilib- 
rium [Fig. [He)]. 

Here we considered two atomic configurations for 
asymmetrically strained bilayer graphene. Because pris- 
tine bilayer graphene has Bernal stacking between top 
and bottom layers [Fig.QJb)], we apply strains along the 
two highly symmetric (i.e., zigzag and armchair) crystal- 
lographic directions of the hexagonal lattice of the bot- 
tom layer. First we fix the stacking order along the short 
edge of the unit cell, and then stretch the bottom layer 
along the normal to the edge. When the bottom layer 
is under strain along the zigzag-chain direction, we call 
this configuration the asymmetrically Z-strained bilayer 
graphene (AZS bilayer) [Fig. QJc)]. When strain is ap- 
plied along the armchair-chain direction in the bottom 
layer, we call it the asymmetrically A-strained bilayer 
graphene (AAS bilayer) [Fig. GJcl)]. We note that, in 
single-layer graphene, strains along those two directions 
induce different variations in its electronic structures and 
the effect of strain along arbitrary directions can be rep- 
resented as a linear combination of them (j. 23]. 

In AZS bilayer graphene, an energy gap opens immedi- 
ately when the bottom layer is stretched along the zigzag 
chain direction; the position where the energy gap opens 
in the first BZ shifts depends on the magnitude of strain 
and is not located at the highly symmetric points of the 
first BZ [Fig. [2]. The fundamental energy gap is indirect 
at strains up to 9%, but direct at greater strains. In AAS 
bilayer graphene, no energy gap opening occurs up to a 
strain of about 14% (not shown here). This shows that 
electronic band gap of asymmetrically strained bilayer 
graphene strongly depends on the direction of strain. 

The energy gap of AZS bilayer graphene is created by 
strain-induced perpendicular electric fields that break the 
symmetry in the onsite energy of the top and bottom 
layers [3, EH . The transverse electric fields across the 
two layers without any external electric source are pos- 
sible because the increase of the work function of the 
strained bottom graphene generates a net charge trans- 
fer from the top layer to the bottom layer [Fig. [3]. The 



FIG. 2: (Color online) The energy contours for (a) the con- 
duction band and (b) the valence band of the AZS bilayer 
with 4% strain on the bottom layer. Only contours for a part 
of the first BZ are drawn, where energy bands exhibit minimal 
and maximal values. Here, TX is one tenth of TX' in the first 
BZ shown in Fig. [TJf ) . The scale bar for the contour is given 
with respect to the Fermi energy (Ef) in units of eV. (c) the 
electronic band structure of the AZS bilayer with 4% strain 
(left panel) along the high symmetric line (rX) and (right 
panel) along the dotted line in (b) that crosses the k-points 
at which the minimum indirect band gap occurs, (d), (e) and 
(f) Energy contours and band structures for the AZS bilayer 
with 6.25% strain using the same convention as in (a), (b) and 
(c). (g), (h) and (i) Energy contours and bands for the AZS 
bilayer with 11.1% strain with the same previous convention. 
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FIG. 3: (Color online) (a) Bilayer graphene in equilibrium 
(left panel) and the effective potential for electrons (right 
panel). Two vertical lines denote the position of bottom 
and top graphene layers, respectively, (b) Asymmetrically 
strained bilayer graphene with the bottom layer uniformly 
stretched along the orange arrows. A uniform electric field 
across the two layers (red arrows) is generated due to work 
function difference (left panel). The effective potential for 
electrons for AZS bilayer with 14.3% strain on the bot- 
tom layer (right) and corresponding electric field (red arrow) 
across the two layers are shown, (c) The work function dif- 
ference (red) and the net potential difference (blue) between 
the top and bottom layers as a function of strain (%) for the 
AZS bilayer. (d) The direct (red) and indirect (blue) energy 
gap of the AZS bilayer as a function of strain. 



electric fields shift and hybridize the conduction and va- 
lence bands, thereby opening energy gaps at the points 
at which band mixing occurs at the Fermi level. The gap 
opening in this system is thus a truly physical manifesta- 
tion of the pseudo-scalar potential induced by mechanical 
deformations of graphene. In our calculations, the work 
function of single layer graphene increases almost linearly 
(0.4 eV by 14% strain) and accordingly, so does the net 
charge transfer between the two layers [Fig. [3fc)]. Due 
to the screening of the transferred charge, the net differ- 
ence between local potentials in the two layers in bilayer 
graphene increases by 0.2 eV for when asymmetric strain 
is 14% [Fig. [3jc)]. To see the change of onsite energies 
of each layer clearly, we calculate the effective potentials 
for electrons in the systems by substracting the total po- 
tential of each single-layer graphene from that of bilayer 
graphene and then taking average in the plane [right pan- 
els in Figs. [3] (a) and (b)]. The induced dipole field across 
the two layers is clearly visible when the lower layer is 
under the strain of 14.3% as shown in the right panel of 
Fig. El (b). 

The indirectness of the energy gap in the AZS bi- 
layer originates from the inequivalent formation of strain- 
induced vector potentials in the two graphene layers. 



In single layer graphene, homogeneous strains shift the 
Dirac cones away from the high symmetric points in the 
first BZ, and this effect is equivalent to introducing a con- 
stant vector potential in the Dirac Hamiltonian [5[l8l.l23|. 
The constant vector potential usually has no effect be- 
cause it can be simply gauged away. In the present case, 
only the bottom layer has such a constant vector poten- 
tial and it cannot be trivial. The displacement of Dirac 
cones in the strained bottom layer alters the points of 
band crossing between the conduction and valence bands 
of the bilayer, which should meet at the K-points of the 
first BZ of bilayer graphene without such asymmetric 
strains (El Hlj|. Our calculations, suggest that asymmet- 
rically strained bilayer graphene is unique in that both 
pseudo-scalar and pseudo-vector potentials induced by 
the strain open the energy gap and determine its charac- 
teristics. 

The competition between various channels of interlayer 
interactions is yet another important factor that deter- 
mines the size of the energy gap. All A AS bilayer s and 
the highly strained AZS bilayer contain large domains 
of carbon hexagons that have almost a perfect overlap 
(AA stacking [30j) between pairs of layers [Fig. DJcl)]. In 
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bilayer graphene with AA stacking, the low energy elec- 
trons behave as those in two interpenetrating Dirac cones 
of single layer graphene and no energy gap forms even 
with strong transverse external electric fields [3l|. Like- 
wise, when the interlayer interactions in the A A stacking 
dominate over all other interactions, the significant dif- 
ferences in the work function between the top and bottom 
layers are less effective in opening energy gaps than the 
cases with the Bernal stacking. For this reason, when 
the strain exceeds 9% in the AZS bilayer, the energy gap 
starts to decrease even though the net potential difference 
continues to increase. Moreover, when strain is applied 
along the armchair direction in the bottom layer of the 
A AS bilayer, the interlayer interaction in the A A stack- 
ing domain is predominant over all other channels and 
the energy gap does not open at all [Fig. DJd)]. 

We find that as the strain in the bottom layer in- 
creases, the energy gap of the AZS bilayer first increases 
linearly, then decreases with a transition from an indi- 
rect gap to a direct gap. At a strain of 5%, both the 
indirect energy gap (0.076 eV) and the direct energy 
gap (0.14 eV) are maximal [Fig. [3(d)] . This variation 
of energy gaps depending on the magnitude of asym- 
metric strains in bilayer graphene is quite in contrast 
to that observed in simple Bernal-type stacked bilayer 
graphene under external electric fields pj], EH, EH, l2o| . 
The latter shows a monotonic increase and then satu- 
ration at an energy gap of ~ 0.25 eV due to screening 
of induced charges 11, 3, 1q| as external fields increase. 
The unique interplay between atomic stacking and strain- 
induced pseudo-scalar potential leads to the maximal en- 
ergy gap at a certain magnitude and direction of strain. 
Hence, even under very gentle strains, the rectifying cur- 
rent due to the strain-induced energy gap should be mea- 
sured in the single-layer/bilayer graphene junctions con- 
sidered here. The atomic orientation of bilayer graphene 
and the size of applied strains could also be inferred from 
such measurements. 

In summary, we show that one can control the en- 
ergy gap of bilayer graphene only by the strength and 
direction of homogeneous strain within the experimen- 
tally accessible range [27|, |32|. We show that, if homo- 
geneous strains with different strengths are applied to 
each layer of bilayer graphene, transverse electric fields 
across the two layers can be generated without any ex- 
ternal electronic sources, thereby opening an energy gap 
and enabling the system to behave as an electromechan- 
ical switch with tunable energy band gaps. This finding 
is a materialization of pseudo-vector and pseudo-scalar 
potentials induced by external mechanical perturbations 
in the relativistic massless Dirac Hamiltonians. 
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